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Abstract 

Triangular mass matrices for neutrinos and their charged partners contain full infor- 
mation on neutrino mixing in a most concise form. Although the scheme is general and 
model independent, triangular matrices are typical for reducible but indecomposable 
representations of graded Lie algebras which, in turn, are characteristic for the stan- 
dard model in noncommutative geometry. The mixing matrix responsible for neutrino 
oscillations is worked out analytically for two and three lepton families. The example 
of two families fixes the mixing angle to just about what is required by the Mikheyev- 
Smirnov-Wolfenstein resonance oscillation of solar neutrinos. In the case of three families 
we classify all physically plausible choices for the neutrino mass matrix and derive in- 
teresting bounds on some of the moduli of the mixing matrix. 
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1. While for many years neither direct searches nor terrestrial oscillation experiments have 
given any positive evidence for nonvanishing neutrino masses, the intense experimental in- 
vestigation of solar neutrino fluxes significantly changed this situation. Indeed, the four solar 
neutrino experiments which test different parts of the energy spectrum and all of which 
find fluxes reduced as compared to the ones expected on the basis of standard solar model 
calculations jlj and non-oscillating electron neutrinos, seem to lend support to the hypoth- 
esis of resonance enhanced oscillation [|[ || within the sun's interior. The flux observed by 
Davies et al. using inverse /3-decay on 37 C1 which is about one third of the expected flux and 
which is sensitive to the high-energy neutrinos from boron, the flux seen in the two exper- 
iments with 71 Ga, GALLEX and SAGE, which is about 60 percent of the expected, and the 
kamiokande results on ^-electron scattering, seem compatible with the Mikeyev-Smirnov- 
Wolfcnstein mechanism calculated from standard solar parameters and two-generation mixing 
with squared mass difference and mixing angle (vacuum values) in the range Q] 

Am 2 w 6 • f (T 6 eV 2 , sin 2 6 w 0.006 . (f ) 

In the meantime further experimental hints at nonvanishing neutrino masses is provided 
by observation of atmospheric neutrinos [^J, and by an accelerator experiment at LAMPF 
. We note that empirical analyses of all five experiments indicate that they cannot all be 
fit with three families of leptons Q . 

Conventional gauge theories of electroweak interactions are based on structure groups 
which are compact Lie groups. Fermions are classified in finite-dimensional, unitary (in gen- 
eral reducible) representations pi of the structure group. As such representations are always 
fully reducible the gauge bosons of electroweak interactions, a priori, can mediate only be- 
tween members of the same, given fermionic family. Any transition between states pertaining 
to two different families must be due to a new element in the theory. 

We adopt the usual attitude in assuming that the mixing of fermionic states with equal 
charges but different flavours is solely due to the mismatch between weak interaction states 
and mass eigenstates. This means that the mixing matrix relevant for the charged current 
(CC) interactions of the standard model is calculable from the mass matrices. Making use 
of the fact that CC interactions involve only left-chiral fields we recently showed that the 
full information contained in any nonsingular mass matrix can be described by means of 
a (lower) triangular matrix The triangular form which is obtained by (unobservable) 
unitary transformations of right-chiral fields, is the most economic and concise description of 
family mixing. As shown in ref. it corresponds effectively to classifying the lepton families 
in reducible but indecomposable representations of the type 

Pi 5 Pi or pi 5 pi 3> pi, (2) 

for two and three families, respectively, where the "semi-sum" means that the right-hand 
representation space is an invariant subspace while the left one(s) is (are) not. 

Although our analysis is completely general, we note that representations of the type 
are not unexpected. Indeed, extensions of gauge theories as obtained in the framework 
of noncommutative geometry add more structure to them and, in particular, constrain the 
fermionic sector in a novel and physically relevant way. In the approach proposed by Connes 
and Lott |j) it is postulated that the Dirac- Yukawa operator contain the physical quark and 
lepton mass matrices, including the quark mixing pattern. The bosonic sector as derived from 



that operator exhibits spontaneous symmetry breaking with the distinctive feature that the 
Higgs potential would vanish if generations did not mix or were degenerate. 

The construction proposed by the Mainz-Marseille group |l(], [llj, on the other hand, 
starts from the underlying graded differential algebra and yields the bosonic sector, including 
spontaneous symmetry breaking, without resorting to the fcrmion content of the theory. 
The Dirac- Yukawa operator is a derived object whose detailed form, unlike the Connes-Lott 
construction, is not an input. The graded gauge potential of this model comprises the gauge 
bosons as well as the Higgs field and is seen to fall into the adjoint representation of the 
graded Lie algebra su(2 | 1), viz. 

[y = g,/ = i] — ► (j = i,y = o)e(/ = o,F = o)e 
(j=i,y = i)ffi(i = ^ = -i). 

The right-hand side gives its decomposition with respect to Lie(SU(2) x U(l)) and de- 
scribes the four gauge bosons, the Higgs doublet and its antidoublet, respectively. As ex- 
plained in earlier publications jy], it is natural to use this algebra as a classifying algebra 
for quarks and leptons as well, its action on the fields of the standard model being under- 
stood in the sense of what we called a weak symmetry in fl2|| . A most remarkable feature 
of su(2 1 1) is the existence of reducible but indecomposable representations which are semi- 
direct sums of two or three copies of the fundamental representations that correctly describe 
the SU(2) x U(l) quantum numbers of quarks and leptons, respectively 0. Such represen- 
tations, in generation space, have a block triangular form and, therefore, provide a natural 
place for generational state mixing. In |l5| we explored this possibility for quark mass ma- 
trices and CKM mixing. Although based on very different assumptions we showed that this 
scheme leads to the mixing pattern favoured by phenomenology, with the bonus of providing 
analytical formulae for the mixing matrix elements in terms of quark masses and a small 
number of free parameters. 

In this letter we study the case of nontrivial neutrino mass matrices, in the general frame- 
work of triangular mass matrices, and their consequences for neutrino oscillations including 
the solar neutrino fluxes. With a pausible physical assumption on the origin of state mix- 
ing, the case of two families contains no free parameters and yields a mixing pattern in fair 
agreement with the resonance oscillation needed to explain the solar neutrino puzzle. The 
case of three families contains more freedom and can be analyzed in terms of complete sets 
of neutrino oscillation experiments. We then illustrate this general framework by the Mainz- 
Marseille model and show that our physics assumptions are quite natural within this model 
based on noncommutative geometry. 

2. When written out in terms of individual neutral and charged fields of definite chirality the 
mass terms are 

^ ass = lti + « } 4 fc) ) + h.c. . (3) 

i=l fc=i 

The matrices T are lower triangular matrices given by 

/ a n \ / dix \ 

= 02i a 22 r<'> = d 2 i d 22 (4) 
\ «3i «32 a 33 J \ d 3 i d 32 d 33 J 

for the neutrinos and the charged leptons, respectively. 



Without loss of generality, weak interaction and mass eigenstates can be chosen such that 
all entries in the triangle matrices (Q) except 031 and ^31 are real jl3||. Changing notations 
slightly we write the latter as 

a 3 ie and d 31 e , 

so that from now on all and all dik are real. The characteristic polynomial for the squared 
mass matrix (T^T^ t) yields the following relations for the parameters 

2 2 2 2 2 2 /r-\ 

m 1 m 2 m 3 = a n a 22 a 33 , (5) 

22, 22, 22 2 2, 2 2, 2 2, 

m 1 m 2 + m 2 m 3 + m 3 m 1 = o, lx a 22 + a 22 a 33 + a 33 a xx + 

a ll a 32 + a 22 a 31 + a 33 a 21 + 

a 2i a 32 _ 2a 22 a 2 ia 31 a 32 cos&"> , (6) 
m\ + m\ + mj = ofj. + af 2 + af 3 + ofi + 031 + 032 ■ (7) 
Similarly, the characteristic polynomial for (T^T^^ yields 

m 2 e mlm 2 T = d 2 n dl 2 dj 3 , (8) 

m 2 e ml + m\m\ + m 2 T m 2 e = d 2 n dl 2 + dj 2 dj 3 + dj^ + 

diid 32 + d 22 d 31 + d 33 d 21 + 

d 2 21 d 2 32 - 2d 22 d 21 d 31 d 32 cos $ (/) , (9) 
m 2 e +ml + m\ = d\ x + d 2 22 + d 2 33 + d\ x + d 2 31 + d\ 2 . (10) 

We now make the following, physical assumption: the mass differences between different 
families are due to electroweak interactions only. This means that without these interactions 
the masses of the electron, muon and tauon would be degenerate, as would the masses of v e , 

and v T . When the interaction is present, its effect within each generation will be the same, 
the observed splittings of the masses will be due to the off-diagonal, inter-family blocks in the 
mass terms (|^), (^). For neutrinos it seems natural to assume further that their primordial 
masses vanish. In terms of the parameters in eq. (^) this means 

an = a 22 = a 33 = , du = d 22 = d 33 = (m e m^m T ) 1/:i . (11) 

As an example consider the case of two families for which the matrices (^) simplify to 

T <» = ( an \ r ( = / d n \ ^ 
V a 2i a 2 2 / \ d 2 \ d 2 2 ) 



These assumptions imply the following choice of the parameters |13 

(a) for the neutrino sector 

an =0 a 22 = a 2 i = ±m 2 , (13) 

which lead to the mass eigenvalues m\ = 0, m 2 ^ 0, and 

(b) for the charged lepton sector 

dn = d22 = ^/m e m^ d 21 = ±(m fl - m e ) . (14) 



The mixing matrix is obtained in close analogy to the case of quark mixing. Let V {v) 
and be the unitary matrices which diagonalize the hermitean products T( V 'T' V " and 
T«T«t, respectively, i.e. V (TT^)V^ = diag ( m 2 ,m 2 ,m 2 ) for {y) and (7). The mixing 
matrix that determines possible oscillations is then 

Vmtx = V^V^K (15) 

In the case of two families it depends on one real angle only whose cosine is 



y/mim e + ^mym^ I WT^ 

cos = = , / — . (16) 

y/(mi + m 2 )(m e + m M ) \jm e +m^ 

From this we obtain the result 

sin 2 (9 = — = 0.0048. (17) 

m e + 

Already this first estimate based on two families of leptons is in fair agreement with the 
value (Q) needed for solar resonance oscillation to occur. The mass eigenvalue of the second 
neutrino remains undetermined from our analysis. The empirical result (jl]) requires it to be 
of the order of 10 -3 eV. 

3. Turning now to the more realistic case of three families we first discuss the neutrino sector. 
Making use of the assumption ( pd| ) eq. (^) implies that at least one neutrino remains massless, 
say, mi = 0, while eqs. (||) - (Q) simplify to 

2222 2,22,2,2 / 1 \ 

m 2 m 3 = a 21 a 32 , m 2 + m 3 = a 21 + a 31 + a 32 . (18) 

Without knowledge of the masses ni 2 and m 3 from some other source there are still various 
possibilities of choosing 021, 031, and a 32 (keeping in mind, of course, that they are not all 
independent). Among these there are two which are particularly interesting. These are 

(A) 031 = 0, where, say, m 2 = a 21 , m 3 = a 32 , 

(B) a 2 i — 0, where, say, m 2 = 0, m 3 — a 31 + a 32 . 

In both cases the matrix is the unit matrix and, by eq. (|l5|), the observable mixing 
matrix is determined solely by the charged leptons, V m ix — V"'t, with V^' expressed in 
terms of mj = m e , m p , m T and of dik as worked out in Setting d\\ = d 22 = d 33 = d we 
find 

/ f(m e )/N e fim^/Nf, f(m T )/N T \ 
V mix = g(m e )/N e g{mJ)/N^ g(m T )/N T , (19) 
V h(m e )/N e /i(m M )/A^ h(m T )/N T J 

where the functions /, <?, and h are given by 

/K) - d 2 d 21 d 32 - dd 3l e~^ il) (d 2 - m 2 ) , (20) 

g{rrii) = m^d^d^e' 1 ^ ' - dd 32 (d 2 - m 2 ) , (21) 

h(mi) = (d 2 - to 2 ) 2 - d 2 21 m 2 . (22) 

The normalization factors in the denominator are 

N e = { [(d 2 - m 2 e ) 2 - m 2 e d 2 21 ] (m 2 - m 2 )(m 2 - m 2 )} 1/2 , (23) 



Nfj, and N T being obtained from this by cyclic permutation of {m e , to m , m T }. Note that due 
to the equations (||) - ( |To| ) there are only two free parameters, say e?2i and d%2- 

The form ( [i"9| ) of the mixing matrix and the formulae - ^3|) are already quite restrictive, 
independently of the choice of e?2i and of d 32 - This is seen as follows. Define the functions 
Fi(z), i = 1,2,3, by 



(m^ — m 2 ){m 2 l — m 2 ) 

with F2 and -F3 obtained from by cyclic permutation of {m e ,m^,m T }. Denoting the 



entries of V m i X , eq. (19), by Vik, we see from eqs. fl22|) and (23) that the moduli of the 



elements in the third line are given by 

\V X \ = Fi(z = dai) , i = 1,2,3. (25) 

Note that because of the ordering of the lepton masses the denominator in F2 is negative 
and, hence, its numerator must also be negative. From these equations one concludes 

K " rf2 |/"V < M21I < (m 2 T - d 2 )/m T . (26) 

When expressed in terms of the r-mass, this means that \d2i\/'m T must lie between 0.0483 
and 0.9993. These boundary values may in turn be inserted into the moduli to obtain limits 
on them, viz. 

0.01006 < IV31I < 0.01116, < IV32I < 0.99995, < |V 33 | < 0.99994. (27) 

Thus, while IV32I and | V33 1 remain practically unrestricted, | V31 1 must always be very small. 

One can also work out limits on the moduli of the elements in the first line of (|l9|) which 
turn out to be even more interesting. In calculating |/(mi)| 2 from eq. (p0[), eliminating d%\ 
and cos $ by means of eqs. (g) and ©, and using d 6 = m^rn^m 2 ., one obtains 

\Vu\ = ( — J Fi(z = d 32 ) , (mi = m e ,m 2 = m»,rri3 = m T ) . (28) 
By the same reasoning as above one sees that ^32 1 must also lie in the interval (p6|). The 



expressions fl28| ) are essentially the same as (g5|), except for the factors (d/m e ) = 89.59, 
(d/rrifj,) — 0.4333, (d/m T ) — 0.02576. Because these factors are so different from one another, 
the allowed intervals ( p7j) for the third line, when applied to the first line, become 

0.9013 < |Vii| < 0.9997, < |Vi 2 | < 0.4333, < \V 13 \ < 0.02576 . (29) 

Thus, while the narrow interval for the first element is shifted to somewhere close to 1, the 
very loose constraints on the second and third element of the third line are turned into rather 
tight limits for the corresponding elements of the first line. 

These limits as obtained on the assumption that neutrinos have no primordial masses and 
yw — 1 (so that V m i X stems from the charged leptons only), could already be tested by 
experiment, either confirming compatibility with this simple scheme or disproving it. Clearly, 
a more detailed analysis will have to include a thorough numerical analysis of acceptable 
values for c?2i and ^32^ and/or other choices of the neutrino mass matrix. 

We note in passing that 1^21 and 1^32 in fact are not completely independent. Using unitarity of V m i x , 
one concludes from the moduli of the elements in the second line that the combination + (d 2 /m^)d| 2 
is bounded from above by 0.9992m 2 , while the combination + (d 2 /m^. )d| 2 is bounded from below by 
0.00287m 2 .. These bounds are so mild that they will not play a decisive role. 



4. As an illustration of this general framework we briefly describe the mass matrices as 
obtained in the framework of the Mainz-Marseille model (cf. introduction). In this model the 
left-chiral doublet and the right-chiral singlet of a given lepton family fit into the fundamental 
representation 

[Y = -l,I = ^]^(I=±,Y = -l)®(I = 0,Y=-2), (30) 

where the right-hand side gives the decomposition in terms of Lie (SU(2) x U(l)). A right- 
handed neutrino would be classified according to the trivial representation [Yq = 0, Iq =0] = 
(I = 0, Y = 0). As discussed in (hJ the two representations can be combined to a reducible 
but indecomposable representation 

Pl := [Y = -1, Jo = 33 [Y = 0, h = 0] , (31) 

which describes a charged Dirac lepton I and its neutrino vf} , along with a right neutrino 
which does not couple directly to the gauge bosons but does couple to the leptonic multiplct 
(Il,Ir, vl) via the Higgs doublet. The next simplest, and physically interesting possibility 
is to classify two or three lepton families in the reducible indecomposable representations 
which are obtained by taking semi-direct sums of the representation pi, eq. (|3l|), according 
to the pattern shown in eq. (|J). While (|3l"l), in a natural way, allows to introduce a (Dirac) 
mass of the neutrino isi, without introducing right-handed interactions, the representation 
(||) provides in addition family mixing both in the charged and in the neutral lepton sector. 
Indeed, in a representation with three families the even and odd generators of su(2 | 1) have 
the form, respectively, 

f, = 



(Ti) 










(Ti) 
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(Ti) 


31 
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(fli 


) n 








ov 


) 21 









)31 


(fii) 32 


(fii) 33 



tk = (Qi) 21 (tt 2 ) 22 . (32) 

V (a) 31 (^) 32 (^) 33 / 

The entries of these block triangular matrices are 4 x 4-matrices. In the basis 

T , , / ,. s ,.,\T 



* = 



(^W,^),^3)) ,with^> = (i« ^4Vj?) , (33) 



and using Pauli matrices Ti and projection operators P± — (1 ± Ta)/2, the odd generators 
must have the form 

(n + r = ( J p aik »- ) (^y k = ( J p _ ) (34) 






CLikP- 


UikP+ 








bikT+ 








<"->"- 1 ^ + "'o + j «)"=( 7 : T _ t). <»> 

the bold face denoting the 2x2 zero matrix. As we showed in |ll| it is no restriction of 
generality to choose the even generators block diagonal, i.e. 

{ht = = (YY k , for (ik) = (21), (31), (32) . (36) 



In the basis chosen the diagonal blocks are 

■5,5,0,0), V2 ( i ± r = 1 



(7 3 r=diag(-i 1, 0,0), V2(7±)»= ( T Z ° ) , (F)« = diag (-1,-1, -2,0). 



The interaction of the lepton state ( J33| ) with the gauge bosons of SU(2) x U(l) and the 
Higgs field is encoded in the superconnection of the model, 

1 - m 1 



A = i lal ■ W + -bYW {8) + - |fi'_^°) + f2' + <£ (+) + ^+^ (0) + ^_0(+)J j , (37) 

with a, b real, dimensionless parameters, /j. the mass scale of the theory. In particular, and as 
explained in ref . jl3| , the mass terms which follow from eq. (^7|) are 

1 n 



r 



(vn+v) + +h. c . , (38) 



where the indices indicate that neutrinos only get contributions from il + , the charged partners 
get contributions only from Thus, in this model C mass has indeed the general form of 
eqs. (§) and (g) above. 

The graded commutators of the even and odd generators of the algebra yield relations 
between the complex parameters a^ . . . 7^ . These reaoQ 

hk = a ik, 0ik = -a»fc) Cik = -dik, lik =Sik, i, k = 1,2,3 , (39) 

a ii 6 i i = 0, z,fc = l,2,3, (40) 

duau = -l, i,fc = l,2,3, (41) 

3 

021^11 + 022^21 = , a 32 (522 + »33^32 = , a 3j£jl = , (42) 

i=l 

3 

£21 an + £22021 = , <5 32 a 22 + <5 33 a 32 = , £3^1=0, (43) 



i=i 

3 



rf2i«ii + ^22^21 = , d 32 a 2 2 + d33«32 = , ^ duotn = . (44) 



i=l 



We note, in particular, that the right-hand sides of eqs. ( |40[ ) and (41) are the electric 
charges of the neutrinos and their charged partners, respectively. The former are important 
in discussing the remaining freedom in choosing the parameters a^ and give an interesting 
interpretation of the cases (A) and (B) discussed above. Indeed, the assumption an = 
complies with the condition ( |40| ) which says that either an or 8n have to vanish because 
neutrinos have no electric charge. Eqs. (^) and ([l3]) then reduce to, respectively, 

a2i£n = , a 32 <522 = 0, a 3 i<5n + a 32 <5 2 i = , (45) 

<>22a2i=0, (533032 = 0, <5 32 a 2 i + <5 33 a 3 i = . (46) 

Barring the trivial case where all are zero, eqs. (^5| - |46|) admit the following choices 
(i) &21 7^ 0, a 3 2 7^ 0: In this case all Sij except <5 3 i must vanish. Two neutrinos are massive, 
one remains massless. Depending on whether a 3 i = the neutrino mass matrix is diagonal 
and the mixing matrix ( |l5| ) is determined by the mass matrix of charged leptons only. 

2 Note that the analogs ojljeqs. (|43[) with 6^. replaced by a^., and a mn by d mn , exist but, by (^l|), are 
linearly dependent on eqs. (|44|). 



(ii) (Z21 ^ 0, 032 = 0: In this case 5n and ($22 must vanish, 821 and J31 remain undetermined, 
while o~32 and £33 are related by 

x ° 31 x 
032 = °33 ■ 

021 

Two neutrinos are massive and mix, as long as 0,31 ^ 0. If 031 = one more neutrino becomes 
massless and, obviously, there is no mixing. 

(iii) a 2i = 0, a 32 ^ 0: Here <5 22 = S33 = 0, 

x 031 x 
021 = on , 

032 

while S32 and S31 remain undetermined. Two neutrinos are massless, the mass matrix is 
diagonal. The mixing matrix ( fl5| ) is determined by the charged sector. 

(iv) CI21 = 032 = 0, but 031 ^ 0: Here Sn = 633 = 0, all others remain undetermined. Two 
neutrinos are massless, the mass matrix is diagonal. 

Note that cases (ii) - (iv) are not simply limiting cases of case (i) because in the latter 
all Sik but S31 must vanish whereas in cases (ii) to (iv) some of these parameters may be 
different from zero. As we noted in sec. 3 the physically interesting cases are (A) and (B) 
which are contained in (i) and (iii). 

5. As CC weak interactions involve left-chiral fields only, triangular mass matrices of the type 
(^) for neutrinos and their charged partners provide the most general, yet concise, description 
of leptonic state mixing. It seems plausible that in the absence of electroweak interactions, 
leptons are mass-degenerate in each charge sector and that neutrinos possess no primordial 
masses. These assumptions are equivalent to the choice ( |TT| ) for the diagonal entries of 
and T®. The characteristic polynomial for the mass matrix (T^T^ t) allows to classify the 
possible choices for the remaining parameters 021, 031, and 032- Among these the choices (A) 
and (B) are particularly interesting because in these cases the neutrino mass matrix is diag- 
onal and, therefore, the mixing matrix is determined by the mass matrix of charged leptons 
only. If the atmospheric neutrino experiment or the short-baseline oscillation experiment is 
confirmed, then case (B) is excluded because one needs (at least) two nonvanishing mass dif- 
ferences. Note that even though, by our assumption of vanishing primordial neutrino masses 
(|ll|), one neutrino stays massless there is nevertheless mixing between all three families. This 
is due to the fact that the electron has a finite mass and, in the charged lepton's mass matrix, 
does not decouple from the other two families. 

The example of two families where all parameters are fixed by the masses, shows that mix- 
ing is compatible with the specific resonance oscillation needed to explain the solar neutrino 
fluxes, cf. eq. (|l7|). With three families there is more, but not much more, flexibility. In the 
case where the mass matrix of neutrinos is such that = 1, i.e. where V m i X depends on 
the charged lepton sector only, there are two free parameters bounded by the interval ( ^6| ) 
(and some mild compatibility condition among them). However, the scheme is so tight that 
even without fixing these parameters, the moduli of the elements in the third and first lines 
of V m i X must lie in the intervals (|27]) and (^), respectively. 

In a more general situation where is not proportional to the unit matrix l the mixing 
matrix ( |l5| ) will also depend on the neutrino masses. We hope to return to a more detailed 
numerical analysis in a future publication. 

Finally, we illustrate this general analysis by a model that is an extension of the standard 
model within noncommutative geometry. This model sheds new light on the classification of 



parameter choices by relating them to graded commutators of generators of the underlying 
graded Lie algebra, and, thereby, establishes a link to lepton quantum numbers. 
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